Incorporating outlier information into diffusion-weighted MRI modeling for robust microstructural imaging and structural brain connectivity analyses

A B S T R A C T The white matter structures of the human brain can be represented using diffusion-weighted MRI tractography. Unfortunately, tractography is prone to find false-positive streamlines causing a severe decline in its specificity and limiting its feasibility in accurate structural brain connectivity analyses. Filtering algorithms have been pro-posed to reduce the number of invalid streamlines but the currently available filtering algorithms are not suitable to process data that contains motion artefacts which are typical in clinical research. We augmented the Con-vex Optimization Modelling for Microstructure Informed Tractography (COMMIT) algorithm to adjust for these signals drop-out motion artefacts. We demonstrate with comprehensive Monte-Carlo whole brain simulations and in vivo infant data that our robust algorithm is capable of properly filtering tractography reconstructions despite these artefacts. We evaluated the results using parametric and non-parametric statistics and our results demonstrate that if not accounted for, motion artefacts can have severe adverse effects in human brain structural connectivity analyses as well as in microstructural property mappings. In conclusion, the usage of robust filtering methods to mitigate motion related errors in tractogram filtering is highly beneficial, especially in clinical stud-ies with uncooperative patient groups such as infants. With our presented robust augmentation and open-source implementation, robust tractogram filtering is readily available.Peer reviewe

On the Viability of Diffusion MRI-Based Microstructural Biomarkers in Ischemic Stroke

Recent tract-based analyses provided evidence for the exploitability of 3D-SHORE microstructural descriptors derived from diffusion MRI (dMRI) in revealing white matter (WM) plasticity. In this work, we focused on the main open issues left: (1) the comparative analysis with respect to classical tensor-derived indices, i.e., Fractional Anisotropy (FA) and Mean Diffusivity (MD); and (2) the ability to detect plasticity processes in gray matter (GM). Although signal modeling in GM is still largely unexplored, we investigated their sensibility to stroke-induced microstructural modifications occurring in the contralateral hemisphere. A more complete picture could provide hints for investigating the interplay of GM and WM modulations. Ten stroke patients and ten age/gender-matched healthy controls were enrolled in the study and underwent diffusion spectrum imaging (DSI). Acquisitions at three and two time points (tp) were performed on patients and controls, respectively. For all subjects and acquisitions, FA and MD were computed along with 3D-SHORE-based indices [Generalized Fractional Anisotropy (GFA), Propagator Anisotropy (PA), Return To the Axis Probability (RTAP), Return To the Plane Probability (RTPP), and Mean Square Displacement (MSD)]. Tract-based analysis involving the cortical, subcortical and transcallosal motor networks and region-based analysis in GM were successively performed, focusing on the contralateral hemisphere to the stroke. Reproducibility of all the indices on both WM and GM was quantitatively proved on controls. For tract-based, longitudinal group analyses revealed the highest significant differences across the subcortical and transcallosal networks for all the indices. The optimal regression model for predicting the clinical motor outcome at tp3 included GFA, PA, RTPP, and MSD in the subcortical network in combination with the main clinical information at baseline. Region-based analysis in the contralateral GM highlighted the ability of anisotropy indices in discriminating between groups mainly at tp1, while diffusivity indices appeared to be altered at tp2. 3D-SHORE indices proved to be suitable in probing plasticity in both WM and GM, further confirming their viability as a novel family of biomarkers in ischemic stroke in WM and revealing their potential exploitability in GM. Their combination with tensor-derived indices can provide more detailed insights of the different tissue modulations related to stroke pathology

Quantitative magnetic resonance imaging towards clinical application in multiple sclerosis

Imaging; Multiple sclerosis; Quantitative MRIImatges; Esclerosi múltiple; Ressonància magnètica quantitativaImágenes; Esclerosis múltiple; Resonancia magnética cuantitativaQuantitative MRI provides biophysical measures of the microstructural integrity of the CNS, which can be compared across CNS regions, patients, and centres. In patients with multiple sclerosis, quantitative MRI techniques such as relaxometry, myelin imaging, magnetization transfer, diffusion MRI, quantitative susceptibility mapping, and perfusion MRI, complement conventional MRI techniques by providing insight into disease mechanisms. These include: (i) presence and extent of diffuse damage in CNS tissue outside lesions (normal-appearing tissue); (ii) heterogeneity of damage and repair in focal lesions; and (iii) specific damage to CNS tissue components. This review summarizes recent technical advances in quantitative MRI, existing pathological validation of quantitative MRI techniques, and emerging applications of quantitative MRI to patients with multiple sclerosis in both research and clinical settings. The current level of clinical maturity of each quantitative MRI technique, especially regarding its integration into clinical routine, is discussed. We aim to provide a better understanding of how quantitative MRI may help clinical practice by improving stratification of patients with multiple sclerosis, and assessment of disease progression, and evaluation of treatment response.C.G. is supported by the Swiss National Science Foundation (SNSF) grant PP00P3_176984, the Stiftung zur Förderung der gastroenterologischen und allgemeinen klinischen Forschung and the EUROSTAR E! 113682 HORIZON2020. F.B. is supported by the National Institute for Health Research biomedical research center at University College London Hospitals. J.W. is supported by the EU Horizon2020 research and innovation grant (FORCE, 668039). D.S.R. is supported by the Intramural Research Program of National Institute of Neurological Disorders and Stroke, National Institutes of Health. A.T.T. is supported by an Medical Research Council grant (MR/S026088/1). S.R. is supported by the Austrian Science Foundation (FWF) grant I-3001. P.S. is supported by the Intramural Research Program of National Institute of Neurological Disorders and Stroke, National Institutes of Health. H.V. is supported by the Dutch multiple sclerosis Research Foundation, ZonMW and HealthHolland

Novel structural-scale uncertainty measures and error retention curves: application to multiple sclerosis

This paper focuses on the uncertainty estimation for white matter lesions (WML) segmentation in magnetic resonance imaging (MRI). On one side, voxel-scale segmentation errors cause the erroneous delineation of the lesions; on the other side, lesion-scale detection errors lead to wrong lesion counts. Both of these factors are clinically relevant for the assessment of multiple sclerosis patients. This work aims to compare the ability of different voxel- and lesion-scale uncertainty measures to capture errors related to segmentation and lesion detection, respectively. Our main contributions are (i) proposing new measures of lesion-scale uncertainty that do not utilise voxel-scale uncertainties; (ii) extending an error retention curves analysis framework for evaluation of lesion-scale uncertainty measures. Our results obtained on the multi-center testing set of 58 patients demonstrate that the proposed lesion-scale measure achieves the best performance among the analysed measures. All code implementations are provided at https://github.com/NataliiaMolch/MS_WML_uncsComment: 4 pages, 2 figures, 3 tables, ISBI preprin

Differences in white matter reflect atypical developmental trajectory in autism: A Tract-based Spatial Statistics study☆

Autism is a neurodevelopmental disorder in which white matter (WM) maturation is affected. We assessed WM integrity in 16 adolescents and 14 adults with high-functioning autism spectrum disorder (ASD) and in matched neurotypical controls (NT) using diffusion weighted imaging and Tract-based Spatial Statistics. Decreased fractional anisotropy (FA) was observed in adolescents with ASD in tracts involved in emotional face processing, language, and executive functioning, including the inferior fronto-occipital fasciculus and the inferior and superior longitudinal fasciculi. Remarkably, no differences in FA were observed between ASD and NT adults. We evaluated the effect of age on WM development across the entire age range. Positive correlations between FA values and age were observed in the right inferior fronto-occipital fasciculus, the left superior longitudinal fasciculus, the corpus callosum, and the cortical spinal tract of ASD participants, but not in NT participants. Our data underscore the dynamic nature of brain development in ASD, showing the presence of an atypical process of WM maturation, that appears to normalize over time and could be at the basis of behavioral improvements often observed in high-functioning autism

Towards contrast-agnostic soft segmentation of the spinal cord

Spinal cord segmentation is clinically relevant and is notably used to compute spinal cord cross-sectional area (CSA) for the diagnosis and monitoring of cord compression or neurodegenerative diseases such as multiple sclerosis. While several semi and automatic methods exist, one key limitation remains: the segmentation depends on the MRI contrast, resulting in different CSA across contrasts. This is partly due to the varying appearance of the boundary between the spinal cord and the cerebrospinal fluid that depends on the sequence and acquisition parameters. This contrast-sensitive CSA adds variability in multi-center studies where protocols can vary, reducing the sensitivity to detect subtle atrophies. Moreover, existing methods enhance the CSA variability by training one model per contrast, while also producing binary masks that do not account for partial volume effects. In this work, we present a deep learning-based method that produces soft segmentations of the spinal cord. Using the Spine Generic Public Database of healthy participants ($\text{n}=267$; $\text{contrasts}=6$), we first generated participant-wise soft ground truth (GT) by averaging the binary segmentations across all 6 contrasts. These soft GT, along with a regression-based loss function, were then used to train a UNet model for spinal cord segmentation. We evaluated our model against state-of-the-art methods and performed ablation studies involving different GT mask types, loss functions, and contrast-specific models. Our results show that using the soft average segmentations along with a regression loss function reduces CSA variability ($p < 0.05$, Wilcoxon signed-rank test). The proposed spinal cord segmentation model generalizes better than the state-of-the-art contrast-specific methods amongst unseen datasets, vendors, contrasts, and pathologies (compression, lesions), while accounting for partial volume effects.Comment: Submitted to Medical Image Analysi

Learn to Ignore: Domain Adaptation for Multi-Site MRI Analysis

The limited availability of large image datasets, mainly due to data privacy and differences in acquisition protocols or hardware, is a significant issue in the development of accurate and generalizable machine learning methods in medicine. This is especially the case for Magnetic Resonance (MR) images, where different MR scanners introduce a bias that limits the performance of a machine learning model. We present a novel method that learns to ignore the scanner-related features present in MR images, by introducing specific additional constraints on the latent space. We focus on a real-world classification scenario, where only a small dataset provides images of all classes. Our method \textit{Learn to Ignore (L2I)} outperforms state-of-the-art domain adaptation methods on a multi-site MR dataset for a classification task between multiple sclerosis patients and healthy controls

Structural-Based Uncertainty in Deep Learning Across Anatomical Scales: Analysis in White Matter Lesion Segmentation

This paper explores uncertainty quantification (UQ) as an indicator of the trustworthiness of automated deep-learning (DL) tools in the context of white matter lesion (WML) segmentation from magnetic resonance imaging (MRI) scans of multiple sclerosis (MS) patients. Our study focuses on two principal aspects of uncertainty in structured output segmentation tasks. Firstly, we postulate that a good uncertainty measure should indicate predictions likely to be incorrect with high uncertainty values. Second, we investigate the merit of quantifying uncertainty at different anatomical scales (voxel, lesion, or patient). We hypothesize that uncertainty at each scale is related to specific types of errors. Our study aims to confirm this relationship by conducting separate analyses for in-domain and out-of-domain settings. Our primary methodological contributions are (i) the development of novel measures for quantifying uncertainty at lesion and patient scales, derived from structural prediction discrepancies, and (ii) the extension of an error retention curve analysis framework to facilitate the evaluation of UQ performance at both lesion and patient scales. The results from a multi-centric MRI dataset of 172 patients demonstrate that our proposed measures more effectively capture model errors at the lesion and patient scales compared to measures that average voxel-scale uncertainty values. We provide the UQ protocols code at https://github.com/Medical-Image-Analysis-Laboratory/MS_WML_uncs.Comment: Preprint submitted to the journa

Harmonizing Definitions for Progression Independent of Relapse Activity in Multiple Sclerosis: A Systematic Review

IMPORTANCE: Emerging evidence suggests that progression independent of relapse activity (PIRA) is a substantial contributor to long-term disability accumulation in relapsing-remitting multiple sclerosis (RRMS). To date, there is no uniform agreed-upon definition of PIRA, limiting the comparability of published studies. OBJECTIVE: To summarize the current evidence about PIRA based on a systematic review, to discuss the various terminologies used in the context of PIRA, and to propose a harmonized definition for PIRA for use in clinical practice and future trials. EVIDENCE REVIEW: A literature search was conducted using the search terms multiple sclerosis, PIRA, progression independent of relapse activity, silent progression, and progression unrelated to relapses in PubMed, Embase, Cochrane, and Web of Science, published between January 1990 and December 2022. FINDINGS: Of 119 identified single records, 48 eligible studies were analyzed. PIRA was reported to occur in roughly 5% of all patients with RRMS per annum, causing at least 50% of all disability accrual events in typical RRMS. The proportion of PIRA vs relapse-associated worsening increased with age, longer disease duration, and, despite lower absolute event numbers, potent suppression of relapses by highly effective disease-modifying therapy. However, different studies used various definitions of PIRA, rendering the comparability of studies difficult. CONCLUSION AND RELEVANCE: PIRA is the most frequent manifestation of disability accumulation across the full spectrum of traditional MS phenotypes, including clinically isolated syndrome and early RRMS. The harmonized definition suggested here may improve the comparability of results in current and future cohorts and data sets

Estimating axon radius using diffusion-relaxation MRI: calibrating a surface-based relaxation model with histology

Axon radius is a potential biomarker for brain diseases and a crucial tissue microstructure parameter that determines the speed of action potentials. Diffusion MRI (dMRI) allows non-invasive estimation of axon radius, but accurately estimating the radius of axons in the human brain is challenging. Most axons in the brain have a radius below one micrometer, which falls below the sensitivity limit of dMRI signals even when using the most advanced human MRI scanners. Therefore, new MRI methods that are sensitive to small axon radii are needed. In this proof-of-concept investigation, we examine whether a surface-based axonal relaxation process could mediate a relationship between intra-axonal T2 and T1 times and inner axon radius, as measured using postmortem histology. A unique in vivo human diffusion-T1-T2 relaxation dataset was acquired on a 3T MRI scanner with ultra-strong diffusion gradients, using a strong diffusion-weighting (i.e., b = 6,000 s/mm2) and multiple inversion and echo times. A second reduced diffusion-T2 dataset was collected at various echo times to evaluate the model further. The intra-axonal relaxation times were estimated by fitting a diffusion-relaxation model to the orientation-averaged spherical mean signals. Our analysis revealed that the proposed surface-based relaxation model effectively explains the relationship between the estimated relaxation times and the histological axon radius measured in various corpus callosum regions. Using these histological values, we developed a novel calibration approach to predict axon radius in other areas of the corpus callosum. Notably, the predicted radii and those determined from histological measurements were in close agreement